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ABSTRACT 



Context. The mass loss from Wolf-Rayet (WR) stars is of fundamental importance for the final fate of massive stars and their chemical 
yields. Its Z-dependence is discussed in relation to the formation of long-duration Gamma Ray Bursts (GRBs) and the yields from 
early stellar generations. However, the mechanism of formation of WR-type stellar winds is still under debate. 
Aims. We present the first fully self-consistent atmosphere/wind models for late-type WN stars. We investigate the mechanisms 
leading to their strong mass loss, and examine the dependence on stellar parameters, in particular on the metallicity Z. 
Methods. We perform a systematic parameter study of the mass loss from WNL stars, utilizing a new generation of hydrodynamic 
non-LTE model atmospheres. The models include a self-consistent treatment of the wind hydrodynamics, and take Fe-group line- 
blanketing and clumping into account. They thus allow a realistic modelling of the expanding atmospheres of WR stars. The results 
are verified by comparison with observed WNL spectra. 

Results. We identify WNL stars as very massive stars close to the Eddington limit, potentially still in the phase of central H-burning. 
Due to their high L/M ratios, these stars develop optically thick, radiatively driven winds. These winds show qualitatively different 
properties than the thin winds of OB stars. The resultant mass loss depends strongly on Z, but also on the Eddington factor T e , and the 
stellar temperature T+. We combine our results in a parametrized mass loss recipe for WNL stars. 

Conclusions. According to our present model computations, stars close to the Eddington limit tend to form strong WR-type winds, 
even at very low Z. Our models thus predict an efficient mass loss mechanism for low metallicity stars. For extremely metal-poor stars, 
we find that the self-enrichment with primary nitrogen can drive WR-type mass loss. These first WN stars might play an important 
role in the enrichment of the early ISM with freshly produced nitrogen. 
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1. Introduction 
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The strong mass loss from Wolf-Rayet (WR) stars fundamen- 
tally affects the evolution, the final fate, and the chemical yields 
of the most massive stars. Its dependence on metallicity (Z) re- 
ceived renewed interest because of the possible importance for 
early generations of massive stars and the ir chemical yields (see 
iMevnet etai]|2006tlChiappini et al.ll2006l). Moreover, within the 
collap sar model ( Wooslev et alj Il993t iMacFadven & Wooslevl 
I1999 ). fast rotating WR stars are now broadly accepted as the 
progenitors of long-duration Gamma Ray Bursts (GRBs). The 
formation of GRB progenitors with the demanded properties 
could b e explained by lowered WR mass loss rates at lo w metal- 
licities dYoon & Langed2005tlWo"oslev & Hegenl2006[) . 

Despite their importance, the mechanism of formation of 
WR-type stellar winds is still under debate. The high WR lu- 
minosities and surface temperatures suggest that WR winds are 
driven by radiation. The observed wind performance numbers 
77 = MVool(L+lc), on the other hand, appear to challenge this 
hypothesis. Commonly observed values in the range 77 = 1- 
10 imply a mechanical wind momentum much in excess of the 
momentum of the radiation field. If WR winds are driven by 
radiation, photons must be used more than once. This can ei- 
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ther be achieved by multiple line scattering, or by successive 
re-distribution and absorption of thermalized photons. 

The question of the driving mechanism of WR mass loss 
is strongly related to its dependence on metallicity (Z). If WR 
winds are driven by radiation, one would expect a similar Z- 
dependence as for the radiatively driven winds of OB star s . Such 
a relation has recently been claimed bv lVink fedeKoten 1 (120051) 
for the winds of late WR subtypes (WNL and WCL). Because of 
the absence of observations over a broad metallicity range, com- 
putational wind models are needed to address these questions. 
However, as previously mentioned, WR winds are far from being 
fully understood. Consequently only a few models are available, 
partly based on controversial model assumptions. 

The first models for WR winds were constructed by 
iLucv & Abbottl (1 1993b and ISpringmannl d 19941) . using Monte- 
Carlo techniques. These authors could show that multiple line 
scattering can account for the high observed wind performance 
numbers (up to 77*10) if enough line opacities are present. 
However, the models could not pove whether WR winds are 
entirely driven by radiation, because they rely on an adopted 
velocity structure, without solving the wind hydrodynamics. In 
particular the models failed to reproduce the wi nd acceleration 
in deep atmospheric layers. The WNL models by de Koter et all 
(119971) and the Z-dependent models bv lVink& de Koterf d2005h 
are based on a similar modelling technique where the ionization 
conditions in the extended WR atmospheres are inferred from 
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model atmosphere grids. Also these models cannot explain the 
wind driving in deep layers. 

A complete ly diff erent approach was f ollowe d by 
INugis & Lamersl (|2002|) and lLamers & Nugisl d2002l) . By 
means of a critical-point analysis these authors could show that 
the observed WR mass loss rates are in agreement with the 
assumption of optically thick, radiatively driven winds. In this 
approach it is assumed that the sonic point is located in very 
deep atmospheric layers where the diffusion limit is valid (note 
that under this assumpti on the sonic point becomes the critical 
point of the wind flow). INugis & Lamersl found that the sonic 
point must be located in specific temperature regimes where the 
Rosseland mean opa city is increasing due to the well-known 
Fe opacity peaks (cf. Ilglesias & Rogerslll996l) . In this picture, 
late-type WR winds are initiated by the "cool Fe-peak" (for 
sonic point temperatures of T s = 40-70 kK), whereas the 
winds of early subtypes are supported by the "hot Fe-peak" 
(T s > 160kK). 

Our own approach follows a computationally much more ex- 
pensive path, where the full wind is modeled in a self-consistent 
manner. For this purpose a hydrodynamic solution scheme is 
incorporated into the Potsdam Wolf-Rayet (PoWR) atmosphere 
models. Within these models a detailed non-LTE radiative trans- 
fer is performed, which takes Fe-group line-blanketing and wind 
clumping into account. They are thus able to describe the con- 
ditions in the WR wind, including the optically thick part, in 
a realistic manner (see Sect. [3] for a detailed description of the 
PoWR models). Using these models we could show that it is 
in principle po ssible to drive a strong W R wind by radiation 
pressure alone (|Grafener & Haman n 2005). In accordance with 
INugis & Lamersl (120021) . we found that the inner wind regions of 
early-type WC stars are supported by line opacities from the hot 
Fe-peak (ionization stages Ferx-Fexvi), and that the observed 
mass loss rates can only be reached for extremely high critical- 
point temperatures (» 200 kK). Moreover, we could demonstrate 
the importance of clumping in the wind dynamics, an effect 
which has not been taken into account in any previous wind cal- 
culation. 

In the present work we focus on late-type WN stars 
(WNL stars, spectral subtypes WN6-9). Compared to early 
subtypes, these stars have rather low effective temperatures 
(30-60 kK) in a sim ilar range as, e.g., typical O supergiants. 
lHamann et al.l d2006l) recently found that these stars form a 
distinct group within the galactic WN population, with ex- 
ceptionally high luminosities (above ~10 5 9 L Q ). This result 
is in line with previous detections of v ery high luminosities 
for H-rich WNL subtype s in the LMC dde Koter et al.l 119971: 
Crowther & DessarJ 19981) . and the galactic center dNaiarro et all 



20041) . Here we investigate the question of why these stars show 



enhanced, WR-type mass loss, and present a grid of WNL wind 
models for different metallicities. Moreover, we examine the 
possibility of strong WR-type mass loss for the earliest gener- 
ations of metal-poor massi ve stars, as a result o f self-enrichment 
with primary nitrogen (see lMevnet et al.ll2006l) . 

In Sect.|2]we give a short overview of the applied modelling 
technique and the adopted atomic data. In Sect.[3]we present the 
results of our model computations and compare them with obser- 
vations. In Sect. [4] we discuss the specific properties of WR-type 
winds, and relate them to previous mass loss studies. Finally, we 
discuss the most important implications of our results in Sect.[5] 
and summarize the main conclusions in Sect. [6] 



2. Atmosphere models 

The Potsdam Wolf-Rayet (PoWR) model atmosphere code 
is a state-of-the-art code for expanding stellar atmospheres, 
which incorporates the treatment of line-blanketing due to 
millions of Fe-group tra nsitions in non-LTE (for details of the 
nume ri c treatment see [Koesterke et al.l 119921: lHamann et al 
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Koesterke et al.l2002tlGrafener et aU2002l: lHamann & Grafene 
20031) . Within our group this code has be en extensively applied 



to the quantitative analysis of WR stars dGrafener et al 
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Hamann et al.l 12006: Ba rniske et al.l 120061) . Recently the code 
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was extend ed by a self-consistent solu tion of the hydrodynamic 
equations dGrafener & Hamannl 120051) . where the atmosphere 
structure (p(r) and v(r)) is determined consistently with the 
radiative acceleration as obtained from the non-LTE radiation 
transport. 

Here we give an overview of the applied modelling tech- 
nique. In Sect. l2.ll we describe the applied numerical methods, 
in Sect. 12. 21 we review the relevant model parameters, and in 
Sect. 12. 31 we summarize the atomic data utilized for the calcu- 
lations. 



2. 1 . Numerical methods 

The model code computes the radiation field, the atomic level 
populations, the temperature structure, the density structure, and 
the velocity field for a stationary, spherically expanding stellar 
atmosphere. The complete solution comprises four parts which 
are iterated until consistency is obtained: the radiation transport, 
the equations of statistical equilibrium, the energy equation, and 
the hydrodynamic equations. 

The radiation transport for a spherically expanding at- 
mosphere is formulated in the co-moving frame of refer- 
ence (CMF), negle cting aberration and advection terms (see 
iMihalas et al.lll976t) . For a fast solution and a consistent treat- 
ment of electron s cattering we employ th e method of variable 
Eddington factors (lAuer & Mihalaslfl970l) . This means that the 
moment equations are solved to obtain the angle-averaged ra- 
diation field, and the numerically expe nsive ray-by-ray trans- 
fer is only calculated from time to time dKoesterke et al.ll2002l: 
iGrafener et alj|2002l) . The fast numerical solution allows a de- 
tailed treatment of millions of spectral lines on a fine frequency 
grid. 

The atomic populations and the electron density are deter- 
mined from the equations of statistical equilibrium. This sys- 
tem of equations is solved in line with the radiation transport 
by applicati on of the ALI forma lism (accelerated lambda it- 
eration, see lHamannlll985l Il986l) . Complex model atoms of 
He, C, N, O, Si, and the Fe-group are taken into account. For 
the inclusion of millions of iron-group trans itions we take ad - 
vantage of the super-level concept following lAndersonl d 19891) . 
The details of the nu merical implementation are described by 
iGrafenere t al. (2002). Density inhomogeneities are taken into 
acc ount in the limit of small-scal e clumps, or micro-clumping 
(see lHamann & Koesterkelll998al and Sects.|4j2]and|53]in the 
present work). The clumping factor D denotes the density en- 
hancement within a clump, or the inverse of the volume filling 
factor fy. 

The temperature structure is obtained from the assumption of 
radiative equilibrium. In the present work this constraint equa- 
tion is decoupled from the equations of statistical equilibrium. It 
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is solved by a temperature c orrection procedure which is based 
on the Unsold-Lucy method dUnsoldll 1955b iLucvll 1 9641) . and has 
been generaliz ed for application in non-LT E models with spher- 
ical expansion (lHamann & Grafenerll2003l) . 

The density- and velocity structure of the expanding at- 
mosphere is obtained from the hydrodynamic equations for a 
stationary radial flow, acc elerated by radiation pressure (see 
IGrafener & Hamannl 120051) . The radiative acceleration is ob- 
tained by direct integration within the CMF radiation transport 



i 4tt r 03 

Orad = K V H V dv, 

P c Jo 



(1) 



where the Eddington flux H v and the opacity k v are both com- 
puted on a fine frequency grid. 

Due to Doppler shifts, the radiative acceleration on spectral 
lines reacts to changes of the velocity gradient v'. Because of its 
influence on the dynamic properties of the gas flow, this depen- 
dence must be included in the hydrodynamic solution. In our 
computations we perform an additional radiative transfer cal- 
culation in advance of each hydro-iteration, where the veloc- 
ity distribution v(r) is changed. In this way we determine the 
fraction of a ra d that responds linearly to small variations of v'. 
This value corresponds to the effective contribution of optically 
thic k spectral lines (a t hick) to the total radiative acceleration (see, 
e.g. JPuls et all 120001) . Dividing a t hick by the total acceleration 
from spectral lines (flii ne s) gives the force multiplier parameter 
a, which has a similar meaning as t he parameter a in t he stan- 
dard theory of line driven winds by ICastor et alj d 1 9751 CAK). 
Note, however, that our a is an effective value which includes 
complex effects like line overlaps and changes in the wind ion- 
ization, while the CAK a is deduced from atomic line statistics. 



2.2. Model parameters 

Our hydrodynamic model atmospheres are characterized by the 
mass, the luminosity, the radius, and the chemical composition 
of the stellar core. Moreover, the clumping factor must be pre- 
scribed throughout the wind. Given these parameters, the density 
and velocity structure of the atmosphere follow from the hydro- 
dynamic solution. The basic model parameters are: the stellar 
core radius /?* at Rosseland optical depth tr = 20, the stellar 
temperature (related to the luminosity L + and the core ra- 
dius R+ via the Stefan Boltzmann law), the stellar mass M*, the 
chemical composition, and the clumping factor D which is pre- 
scribed as a function of tr. Note that we define tr in this context 
as the 'LTE-continuum' Rosseland mean opacity. 

For the line broadening we assume Doppler profiles with a 
broadening velocity of i>d = lOOkms -1 throughout this work. 
In this way the absorption troughs of P-Cygni type line profiles 
are successfully reproduced. Moreover, we are presently limited 
to relatively high broadening velocities because the frequency 
spacing in our models scales with uu, and the computation time 
thus scales with 1/i>d- 

Because of its effect on the radiative acceleration, the de- 
tailed prescription of the clumping factor D(r) has an influence 
on the obtained results. Unfortunately, this value can only be 
constrained very roughly from the obse rved electron scatter- 
ing wi ngs of strong emission lines (see lHamann & Koesterkel 
Il998al) . For the outer part of the wind we thus adopt a constant 
value of Z) max = 10, which ga ve satisfactory resu l ts in previous 
analy ses of similar stars (e.g.. iHerald et al.ll2.Q01b iMartins et aT] 
120081) . For the radial dependence in the deeper wind layers we 
rely on previous results for WR 111 ( Graf ener & Hamannl2005l) 
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Table 1. Summary of the model atom. Fe-group ions (Fe) are de- 
scribed by a relatively small number of super-levels, each repre- 
senting a large number of true atomic energy levels (sub-levels). 
The relative abundan ces of the Fe-group elements are listed in 
IGrafener et al. (2002), Table 2. 



where a r-dependent formulation reproduces high ionization 
stages very well, which are formed deep inside the atmosphere. 
In this formulation D continuously increases from a value of 
D m i n = 1 at the wind base (for tr > 0.7) to D max = 10 in the 
outer layers (for tr < 0.35). 

The freedom in the prescription of D(r) is not satisfactory. 
However, clumping is clearly detected in WR winds and should 
be included in the models. As long as only optically thin regions 
above the critical point of the equation of motion are affected 
by clumping, its influence on the obtained mass loss rates is rel- 
atively small (see Sect. l4.2l i. Nevertheless, we expect a strong 
influence on the outer wind structure and in particular on the 
terminal wind velocities. 



2.3. Atomic data 

The chemical composition of the model atmospheres is 
given by the mass fractions X He , Xq, X n , Xq, X&, and 
of helium, carbon, nitrogen, oxygen, silicon, and iron- 
group elements. Compared to our pre v ious models for early- 
type WC stars (IGrafener et alj l2002t IStasihska etail l2004t 
IGrafener & Hamannl 120051) we have reduced the model atoms 
significantly. This is possible because of the lower effective tem- 
peratures, and the specific chemical composition of WN stars 
(mainly H, He, and N). In particular, we use rather small atomic 
models for carbon and oxygen, and a reduced number of ioniza- 
tion stages for the iron group (Fei-x). In this way we manage 
to combine the most relevant atomic species in a model atom 
with 410 super levels. This model atom still allows moderate 
computing times of the order of 1-3 days for a hydrodynamic 
model, i.e., it permits parameter studies with a large number of 
grid models. 

The model atoms are summarized in TableQ] The atomic 
data are compiled from the following sources. Oscillator 
strengths for CNO are ta ke n from the O pa city Project 



dHummer & Mihalasl 1 1 988b 



1992b ICunto & Mendozaf 



Seatonl 
11222 



1987 



Seaton et al.l 
ISeaton et al.1 1 1994b 



The Opacity Project Teamlll995l). level e nergies from Kurucz's 
CD-ROM No. 23 dKurucz & Belli 1 19951). and i onizati o n cros s 
sections from iNahar & Pradhanl (fl997h and iNahari dl999l) . 
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Fig. 1. Wind models for galactic WNL stars: mass loss rates for 
different stellar temperatures and Eddington factors T e . The 
corresponding spectral subtypes are indicated tog ether with WR 
numb ers of specific galactic objects (according to Ivan der Huchtl 
1200 ll in brackets), which show a good agreement with the syn- 
thetic line spectra. Note that the models are computed for a 
fixed stellar luminosity of 10 63 L Q . The WN7 model (WR22) 
is slightly offset from the standard grid models because it is cal- 
culated with an enhanced hydrogen abundance (see Sect |3.21 >. 



Fe-group data are ta ken from Kurucz CD-ROM No. 20-22 (see 
lKuruczlll99 11120021) . 



3. Model calculations 

In this section we present a systematic parameter study of the 
mass loss from late-type WN stars. The first question concern- 
ing these objects is -why they become WR stars, i.e., why they 
show enhanced mass loss. We address this question in Sect. 13. II 
where we investigate the dependence of WNL mass loss on stel- 
lar parameters, and compare our results with observed properties 
of the galactic WNL sample. In Sect. |3.2| we perform a more de- 
tailed comparison with WR 22, a spectroscopic binary for which 
the mass of the WR component is roughly known. Finally, we in- 
vestigate the Z-dependence of WNL star mass loss in Sect. |3.3l 
and combine our results in the form of a mass loss recipe in 
Sect.r 



3.1. WNL stars at solar metallicity 

A l arge part of the g alactic WN population has been re-analyzed 
by lHamann et al.l d2006l) . based on a publicly accessible 
grid o f line-blanketed atmosphere models ( Hamann & Grafener] 
|2004 http://www.astro.physik.uni-potsdara.de). In 
this study the WN paramete rs have been si gnificantly revised 
(e.g., compared to lHamann &Koesterkelll998bl) . chiefly due to 
the influence of Fe-group line-blanketing, and the considera- 
tion of 2MASS IR photometry. In particular, the galactic WN 
stars are now found to form two distinct groups which are di- 
vided by their luminosities. The first group consists of a mixture 
of early- to intermediate WN subtypes with luminosities below 
10 5 - 9 L . The second group, formed by the WNL stars (WN 6-9), 
shows luminosities above this value. While the earlier subtypes 
are dominated by pure He-stars, WNL stars tend to show a sig- 
nificant amount of hydrogen at their surface (up to Xu = 0.55). 
They are located to the right of the zero-age main-sequence in a 
temperature range of = 35-55 kK, and thus might still be in 
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Fig. 2. Wind models for galactic WNL stars: wind den s ities c om- 
pared to observed values according to lHamann et al.1 d2006l) . In 
the R t -T+ plane, our models (solid curves) reproduce the upper 
part of the observed WNL sample (crosses), corresponding to 
WNL stars with low wind densities (see text for the definition of 
R ( ). The objects belonging to the enigmatic class of WN 8 spec- 
tral subtypes are indicated in red. 
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Fig. 3. Wind models for galactic WNL stars: ter minal wind ve- 
locitie s compared to observed values taken from lHamann et all 
(120061). and the observed trend with spectral subtype (from 
iNiedzielski & Skorzvnskil [20021) . Our models (solid curve for 
F e = 0.55 and dashed curve for F e = 0.67) reproduce the ob- 
served trend of increasing Uco with increasing T*, but they tend 
to underestimate the high values observed for the hottest objects. 



the phase of central H-burning. In the present work we concen- 
trate on this class of extremely luminous evolved stars, which are 
in many cases among the most luminous objects in young stellar 
clusters. 

Our main goal now is to show why these stars become 
WR stars, or in other words, why they show enhanced mass loss. 
WNL stars have similar effective temperat ures and surface com- 
positions as typical O-supergiants (see e.g. lCrowther et al1l2002l) 
but they show much stronger wind emission lines. To investi- 
gate their mass loss properties we have prepared a series of hy- 
drodyna mic PoWR models w ith typical WNL parameters (re- 
ferring to lHamann et al.ll2006l) . We adopt a high luminosity of 
Li, = 10 63 L Q , effective core temperatures in the range = 33- 
55 kK, a clumping factor of D max = 10, and typical WN surface 
abundances with X H = 0.2, X He = 0.78, and X N = 0.014 (for 
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10 -4.85 


1Q-(4.6...5.4) 


1Q-(4.7...5.5) 


1 Q-(4.2...5.2) 


|Q-(4.2...5.5) 


10 -(5.41) 


km s 


778... 1298 


974 


240... 1903 


388. ..478 


646 ...1396 


620... 1567 


941 




9.6... 20.6 


21.7 


9.2 ...53.7 


9.1 ...37.7 


12.5 ...39.7 


8.02... 30.7 


28.2 


n 


0.28... 1.51 


0.34 


0.04... 1.15 


0.04. ..0.21 


0.20... 1.12 


0.20... 1.72 


0.44 



Table 2. Model parameters for our different (grid) computations. The first column (Grid 1) denotes the models from Sect. l3.1l for 
solar Z, where and F e are varied. In the second column the parameters for the WR22 model from Sect. l3.2l are given. The third 
column (Grid 2) indicates our Z-dependent models from Sect. 13. 31 followed by the computations for extremely metal-poor stars with 
zero Fe abundance (Grid 3). The next two columns belong to Sect. 13. 41 where the dependence on (Grid 4) an d L» (Grid 5) are 
investi gated. The last column (VdK05) indicates a model from Sect. |4.3| for comparison with the WNL model from lVink & de Koterl 
d2005l) for solar Z. Note that for the Z-dependent models z = Z/Z Q indicates the relative mass fraction of metals with respect to 
the solar value. The luminosity L*, the radius the temperature T*, and the mass M* of the stellar core (with the corresponding 
Eddington factor F e ), the assumed maximum wind clumping factor £> max , the surface mass fractions X He , Xc, X H , Xq, X^, and X Fe , 
and the resultant mass loss rates M, transformed radii R t , terminal wind velocities y M , and wind efficiencies r\ - Mv^liL+lc) are 
given. 



more details see Table[2j grid). For our hydrodynamic mod- 
els it is moreover necessary to prescribe the stellar mass. We 
study two values, M* = 67 M and 55 M Q , corresponding to rel- 
atively large Eddington factors F e = 0.55 and 0.67 (note that we 
are giving here T e = ^ e L*/47rcGM* for a fully ionized plasma). 

The obtained mass loss rates are plotted in Fig.Q] We find a 
strong dependence on F e (or equivalently the L/M ratio), and 
on the stellar temperature (M oc Note that both effects 

mark an important change in the properties of WR-type winds, 
compared to the radiatively driven winds of OB stars. As we will 
discuss in more detail in Sect. 14.11 this is a direct consequence 
of the optically thick wind physics, where the location of the 
critical point is coupled to T e and the local electron temperature 
T e . 

The synthetic spectra obtained from our calculations are con- 
sistent with observed WNL spectra. In Fig.[TJwe have indicated 
four specific objects which are in good qualitative agreement 
with our models. They reflect the observed sequence of WNL 
spectral subtypes, starting with WN 6 at 55 kK to WN 9 at 3 1 kK. 
Note that these objects show similar line spectra as our models, 
but different luminosities and surface compositions. Because of 
this, their actual mass loss rates may differ from the values given 
in Fig.Q] A detailed comparison with the WN 7 subtype WR 22 
is presented in the next section. 

For a more comprehensive overvi ew we compare our re sults 
to the stellar parameters obtained by lHamann et al.l (12006b for 
the galactic WNL sample. However, because our computations 
are restricted to a fixed luminosity, we cannot simply compare 
the obtained mass loss rates. We thus employ the transformed 
radius 

Uoo / VDM f /3 
flt=A* t r 7 . (2) 

2500^8-' / 10- 4 M yr-' J 



which is a lumi nosity-independent measure of the wind density. 
As outlined by lHamann & Koesterkd dl998al) . models with the 
same R t are scale invariant, i.e., for a given value of T* and 
a given surface chemistry they show the same line equivalent 
widths. The results are shown in Fig. [2] In the R,-Ti, plane our 
wind models reproduce the upper half of the observed WNL 
sample, i.e., WNL stars with low wind densities. The large scat- 
tering in the observations makes it difficult to draw clear con- 
clusions from this comparison. Note, however, that e.g. constant 
mass loss rates would show a clear trend with R, oc oc 
The fact that the observations as well as our models show no 
clear trend thus supports our finding of increasing M with de- 
creasing r*. The large part of the high wind density objects, 
which are not reproduced by our models, belong to the enig- 
matic WN 8 subtypes. These stars show strong photometric vari- 
ations. Pulsations thus might play a role in their mass loss (see 
the discussion in Sect. l5.4b . Moreover, most of them are located 
far away from stellar clusters, i.e., their luminosities are largely 
unknown. They might thus belong to a different group of objects 
with much lower luminosities. 



In Fig. [3] we compare the obtained terminal wind veloc- 
ities with observations. The galactic WR stars are known to 
show a trend of increasing Vqq for earlier spectral subty pe s 
dHowarth & Schmutzlll992l:lNiedzielski & Skorzvnskil2002l) . In 
Fig. [3] this trend is indicated by the grey shaded area, together 
with o bserved values for single WNL stars from lHamann et all 
(120061) . Our models reproduce approximately the mean of the 
observed values. However, although they show increasing Vao 
with increasing T*, this trend is less pronounced than the ob- 
served one. Particularly around = 50 kK our models tend to 
saturate whereas the observed values still increase. In Sect. 14. 21 
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3.6 3.8 

iog 10 (vA) 

Fig. 5. WR 22: observed flux distribution in absolute units (blue), 
including optical + 2MASS IR photometry (blue boxes), com- 
pared with the model flux corrected for interstellar extinction 
(red). The detailed model parameters are given in Table|2] 



we will show that this issue might be related to the detailed radial 
dependence of the clumping factor D(r). 

3.2. Test case: WR22 

To verify our assumption of high L/M ratios for WNL stars 
we perform a more detailed comparison with WR 22, an eclips- 
ing WR+O binary system in Carina OB 1, which is classi- 
fied as WN 7h. WR 22 is an ideal target for our purposes be- 
cause mass estimates are available for the WN component. 
Moreover, the Ostar component is so faint that at least 90% 
of the observed flux originates from the WR star, i.e., the 
WN spectrum is o nly marginally contaminate d. From the ra- 
dial velocity curve ISchweickhardt et alj ( 19991) derived a m ass 
of M W r sin 3 i = 55.3 + 7.3 M G , whereas iRauw et al.l (1 19961) ob- 
tained 71.7 ± 2.4 M Q , both with M WR /M Q = 2.7. The inconsis- 
tency is presumably due to the high eccentricity of the system, 
and the extremely high signal-to-noise ratio that is need to detect 
the O star absorption lines. 

In Fig. [4] we compare our best fitting synthetic spectrum 
to observations of WR 22. The stellar parameters are the same 
as for our previous gri d models, but with T * = 44.7 kK and 
Xh = 0.4, as derived bv lHamann et al.l (120061) . With our hydro- 
dynamic models we obtain a good spectral fit for an Eddington 
parameter of F e = 0.55, corresponding to a mass of 78.1 M 
(the detailed model parameters are summarized in Table|2]i. The 
resulting model reproduces the observed emission line spectrum 
reasonably well, although the terminal wind velocity appears too 
small and the He i lines are slightly too weak. Nevertheless, the 
fit quality is more than satisfactory for our present purposes. The 
absolute flux distribution of WR 22 (from IR to UV) is closely 



reproduced with a lu minosity of 10 L for an adop ted distance 
modulus of 12.1 mag dLundstrom & Stenholmll 19841). and an in- 
terstellar extinction law according to iFitzpatrickl (1 19991) with 
E B _ V = 0.42 and R = 3.6 (see Fig.g). 

The mass of 78 . 1 M Q , as obtained from our mod els, is appar- 
ently i n good agreement with the mass estimate by iRauw et all 
d 19961) . However, for a proper comparison we need to take into 
account the inclination of the system, its distance, and possiple 
uncertainties in our spectral analysis. 

The inclination of the system is constrained by the detection 
of an ecli pse near periastron , where the O s tar is obscured by the 
WR star (IGosset et alii 19911) . According to ISchweickhardt et all 



(1 19991) the periastron angle of the O star is 88°, and the or- 
bit is highly eccentric with e = 0.598. No eclipse is detected 
near apastron, presumab ly because of the eccentricity of the 
orbit (IGosset et al.l 1 199 ll) . According to our models the radius 
of the stellar disk, i.e., the radius where tr oss = 2/3, is only 
slightly larger than the stellar core radius: R2/3 = 1.056/?* = 
25.0^0. Adopting a sin; = 2.31 10 8 km = 13.3 R 2 p from 
ISchweickhardt et all d 19991) . we obtain a projected distance be- 
tween both stars of 5.35 R2/3 at periastron and 21 .3 R2/3 at apas- 
tron. The eclipse detection at periastron thus constrains sin i to 
values above 0.983 (sin 3 i > 0.950). The non-detection at apas- 
tron limits sin; to values below 0.9982 (sin 3 / < 0.9946), i.e., 
sin 3 i = 0.972 ± 0.022. Note that we have neglected the size of 
the O star in this estimate because a significant part of the O star 
is obscured but it is not clear whether the eclipse is full or partial. 
Applying our limits for sin 3 i we obtain Mwr - 5 6 .9 + 8 .5 M B 
using the orbital data from ISchweickhardt et al l (1 19991) . and 
M W r = 73.8+4.0 M Q based on lRauw et al.l dl996l) . Hence the in- 
clination induces only small additional uncertainties in the mass 
estimate. 

The emission line strengths computed in our models depend 
strongly on the wind density, and thus on the obtained mass loss 
rates. These are chiefly affected by and r e . While T+ can be 
reliably determined from spectral line diagnostics, T e depends 
on the L/M ratio and on the hydrogen abundance X^. Again, 
Xu can be determined from spectral line ratios. Our models 
thus constrain the L/M ratio, with an error of about 10% (see 
Sect. 14. 21 for a discussion of additional systematic errors). The 
obtained value for M* thus scales with L+, which depends on 
the adopted distance and interstellar extinction. The extinction 
is determined within our spectral analysis from the continuum 
shape in the UV. Together with this might introduce an ad- 
ditional error of about 10% in L*. More critical is the adopted 
distance modulus for Cari na OBI. Values in the literature span 
a range between 1 1 .8 mag dDavidson et alj|2001l) and 12.55 mag 
dMassev & Johnsonll 19931) . corresponding to an uncertainty of a 
factor of two in L+. The unknown distance of the system is thus 
the main source of uncertainty in our analysis. 

We have seen that, with a n intermediate distance modu - 
lus of 12.1 mag (according to iLundstrom & Stenholml 1 19841) . 
our model i s in a greement with the mass determination by 
IRauw et al.l (119961) . For the low value of 11.8 mag, WR22 
would become slightly fainter (10 618 L ) and our mass estimate 
would roughly scale down to 59 M , now in agreement with 
ISchweickhardt et all ( fl999h . With 12.55 mag we would obtain 
a very high luminosity of 10 6 5 L which puts WR 22 extremely 
close to the Eddington limit or even above (r e = 0.92 for 74 M , 
and F e = 1 .24 for 55 M ). Note that in view of the additional b-b 
and b-f contributions to the mean opacity also F e = 0.7 would 
imply a super-Eddington luminosity (see Sect. l4.2l i. Such a large 
distance to WR 22 thus seems to be very problematic. 

At the present stage we conclude that WR22 has a very 
high luminosity of the order of 1O 6 3 L , and that it is located 
rather close to the Eddington limit. Chiefly because of the un- 
certain distance, the possible error in the luminosity is large. 
Nevertheless, mass estimates from radial velocity measurements 
support our assumption of large Eddington fa ctors for WNL sub- 
types. With the large distance determined bv lMassev & Johnson! 
(119931) . WR22 would even exceed the Eddington limit. 



3.3. Grid calculations for different metallicities 

The wind acceleration in our models predominantly originates 
from Fe line opacities. We thus expect a strong dependence of 
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Fig. 4. WR22: Comparison of the synthetic spectrum from our hydrodynamic model (dashed line, red) with observations (blue). 
Fluxes are given in units of the model continuum flux. Note that only the optical and IR data are rectified by hand, whereas the 
observed UV fluxes are divided by the model continuum corrected for interstellar absorption and distance. The detailed model 
parameters are given in Table|2] 



the mass loss on metallicity. To investigate this dependence we 
have prepared a grid of hydrodynamic WNL models for a wide 
range of Z. Because we already know about the fundamental 
importance of the Eddington limit for WR-type winds we have 
introduced F e as an additional free parameter. The stellar core 
parameters (L*, T+, /?*), and the surface abundances Xn and Xn e 
are adopted from our WR 22 model in the previous section. M* 
is adjusted to match the desired values of F e , and Z is varied by 
scaling the surface mass fractions Xq, Xn, Xq, Xsi, and Xp e (see 
Table!] Grid 2). 

The results are presented in Fig. [6] The solid (blue) curves 
indicate model series for different initial metallicities, i.e., with 
CNO abundances corresponding to a secondary helium and ni- 
trogen enrichment by the CNO-process. As expected, the models 
show a strong Z-dependence. However, our second free parame- 



ter F e turns out to be at least equally important. In fact it seems 
that the proximity to the Eddington limit is the primary reason 
for the enhanced mass loss of WR stars. Remarkably, the high 
WR-type mass loss rates can even be maintained for very low Z, 
if the stars are close enough to the Eddington limit. 

The dashed (red) curve in Fig. [(^indicates models with Xp e = 
and Xsi = but Xcno = Z (see Table|2j Grid 3). These mod- 
els show to what extent the primary production of nitrogen may 
affect the mass loss at different metallicities. Such an enhance- 
ment is e xpected for fast rotat ing stars at extremely low metallic- 
ities (see iMevnet et alj |2006i). Interestingly, these models show 
a similar behavior to the Fe-rich models. The solar-like CNO 
mass fraction produces a stellar wind that is similar to an Fe- 
rich model with 1/50 Z . This means that the potential mixing 
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Fig. 6. WNL star mass loss in the range 
from Z o /1000 to 2Z Q .: Mass loss rates 
(top) and terminal wind velocities (bottom), 
as obtained from our hydrodynamic grid 
models, are plotted vs. the Eddington fac- 
tor F e (the models are computed for a fixed 
value of = 10 6 3 L Q , the variation of F e 
thus effectively corresponds to a variation 
of the stellar mass). Model parameters are 
given in Table[2] The solid blue lines indi- 
cate model series where T e is varied for a 
given value of Z. Note that strong WR-type 
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if T e approaches unity. Moreover, the mod- 
els show a clear trend of decreasing d m with 
decreasing Z. The dashed red lines indicate 
models with zero Fe, but solar-like CNO 
abundance. They show that strong WR-type 
mass loss is also possible for extremely 
metal-poor stars due to the self-enrichment 
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of primary nitrogen to the surfaces of metal-poor stars consider- 
ably affects the expected mass loss rates. 

In the bottom panel of Fig. [6] we show the terminal wind 
velocities obtained from our models. First, the models clearly 
predict decreasing wind velocities with decreasing Z. This find- 
ing is in agreement with observations although the o bservational 
evidence is not as cle ar as our prediction (see e.g. IConti et all 
Il989t ICrowfherifcOOOl) . Second, the terminal velocities tend to 
stay constant or even increase with M. This behaviour again in- 
dicates important qualitative differences to OB star winds. For 
OB stars the well-established wind momentum-luminosity rela- 
tion predicts the opposite behav iour, i.e., decreasing w ind veloc- 
ities for increasing mass loss dKudritzki et al.| [l999). A closer 
inspection of our models shows that the changes in the velocity 
structure are related to changes in the wind ionization. For cases 
where an increase of M leads to recombination of Fe, the veloc- 
ity increases in the recombination region due to the newly ex- 
posed line opacities. Models with a similar ionization structure, 
on the other hand, show similar wind velocities. Note that our 
assumption of a constant Doppler broadening (uu - lOOkms -1 , 
see Sect. 12. 2t might affect the results for extremely low Z, be- 
cause for these cases the wind velocities are of the same order of 
magnitude as vd- 
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Fig. 7. The dependence of WNL star mass loss on the hydrogen 
abundance X H : for otherwise fixed stellar parameters the mass 
loss depends only weakly on Xu, with M oc Xu 45 . Model com- 
putations for three different values of and fixed values of 
T e = 0.55 and L* = 10 6 3 L o are shown. Note that the stellar 
mass is varied to keep T e constant. 



3.4. A mass loss recipe for WNL stars 

A primary goal of our present work is to provide a recipe that de- 
scribes the mass loss of WNL stars, dependent on the fundamen- 
tal stellar parameters M*, T*, and the surface mass fractions 
Xh, Xue, and Z. In the preceding sections we emphasized that, 
for fixed L*, the mass loss is chiefly determined by the proxim- 
ity to the Eddington limit, i.e., by the Eddington factor F e . We 
thus expect that the main influence of the surface mass fractions 



Xu and Xn e is due to the number of free electrons, which enters 
the Eddington parameter for a fully ionized plasma in the form 
T e = io- 4 813 (l + X H ) In Fig.|7]we show that for fixed T e 

the influence of the light elements is indeed only moderate, with 
log(M) oc -0.45 Xu. As we will discuss in Sect. |4.3l the slight 
decrease of M with increasing hydrogen abundance is due to the 
influence of the He n continuum in the critical layers. 

We present in Fig.[8]a series of models with varied L*. Again, 
the stellar parameters are adopted from WR 22. T e is fixed to a 
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According to Eqs. (0 and (j4]i the WNL mass loss rates would 
go to zero when T e approaches To. Our recipe of course loses 
validity before this point, when the winds become optically thin. 
For our present grid models this happens at a mass loss rate of 
log(M/ M Q yr 1 ) * -5.5 when F e - T » 0.1. Note that this cri- 
terion should roughly depend on the wind density at the sonic 
point. We thus expect the limiting mass loss rate, above which 
our models are valid, to be of the order of log(Mn m / M yr _1 ) * 
-5.5 + 21og(i?*/23.7i? Q ). Note that Eqs. © and © which de- 
scribe our parameters Fo(z) and /3{z) are strictly limited to the 
investigated parameter range of 2 Z Q < Z < 1O~ 3 Z . 



4. The properties of WR-type stellar winds 

In the present section we discuss the properties of WR-type stel- 
lar winds as opposed to the thin winds of OB stars (Sect. 14. it . 
In this context we discuss possible error sources, with emphasis 
on the unknown turbulent velocities and wind clumping factors 
(Sect. l4.2i i. Moreover, we compare our results with other mass 
loss predictions for WR and OB stars (Sect. 14. 3] ). 



value of 0.55, and L* and T+ are varied over the relevant param- 
eter range (see Table|2]i. Notably, the resultant mass loss rates de- 
pend only weakly on L+, with M oc L® A1 . Note that this again de- 
notes an important difference to OB star winds for which a much 
stronger dependence is expected (see the discussion in Sect. K.U . 

We combine our results in the following form. For the de- 
pendence on T+ we use the scaling relation from Sect. |3.1| with 
M oc r^ 3 5 , and for the dependence on Xu and L* we use 
M oc Xu° 45 L° A1 from this section. For the coupled dependence 
on F e and Z from Sect. 13. 31 we use a more complex description 
of the form M oc (F e - Tof, where the Z-dependent parame- 
ters Fo(Z) and /3(Z) are determined by a multi-dimensional fit 
procedure. The resultant mass loss recipe is in good agreement 
with our grid computations in Fig© with an accuracy better than 
0.05 dex except for two data points with a maximum deviation of 
0. 1 dex. Combining all these relations, we arrive at a mass loss 
recipe of the form 

log^wNL/Moyr 1 )^ -3.763 

+j3 ■ log (F e - r ) - 3.5 • (log(TVK) - 4.65) (3) 
+ 0.42 ■ (log(L*/Lo) - 6.3) - 0.45 • (X H - 0.4) , 

or 

log^wNL/Moyr 1 )^ 10.046 

+/Mog(r e -r )-3.51og(7yK) (4) 
+ O.421og(L*/L )-O.45X H , 

with 

yS(Z) = 1.727 + 0.250 ■ log(Z/Z ), (5) 
r (Z) = 0.326 - 0.301 • log(Z/Z Q ) - 0.045 ■ log(Z/Z ) 2 . (6) 

This recipe is only valid for optically thick winds that are driven 
by the cool Fe opacity -peak. It is thus only applicable in a limited 
range of stellar temperatures. should never exceed a value of 
« 70 kK because then the opacity minimum between the hot and 
the cool Fe-peaks is reached. On the cool side, for < 30 kK, 
our models reach a limit where the radiative force exceeds grav- 
ity in the deep hydrostatic layers due to the hot Fe-peak. At this 
point the region of LBV-type instabilities might be reached (see 
the discussion in Sect. |5.4t . 



4. 1 . Optically thick stellar winds 

In Sects. [3~T1 and I3.4I we have seen that our WR wind models 
are qualitatively different from OB-type stellar winds. In par- 
ticular we find a strong dependence on the Eddington factor 
F e , a steep increase of M with decreasing effective temperature 
(M oc r^ 3 5 ), and a weak dependence on the stellar luminosity 
(M oc Our results from Sect. I3.3I even suggest a violation 
of the wind momentum-luminosity relation, whi ch is extremely 
well- e stablished for OB s t ar win ds (see, e.g., iKudritzki et al.l 
I1999I: iMokiem et al.ll2005l l2007bl) . Instead, the velocity struc- 
ture in our models seems to be dominated by ionization effects. 

Within the cla ssical CAK theory for OB star winds 
dCastor et alj 1 1 975b a mass loss relation of the form M oc 
M^r 1 ^" lJ" is expected. As explained in Sect. 12. II the parame- 
ter a characterizes the response of the radiative acceleration to 
the velocity gradient, and typically lies around a as 0.6. M e ff de- 
notes the effective stellar mass, reduced by the contribution of 
the radiative acceleration on free electrons M s g = M*(l - r e ). 
In contrast to our WR models, the expected dependence on the 
stellar mass is thus rather weak, while the luminosity depen- 
dence is comparatively strong. Moreover, OB star winds do not 
have a pronounced wind ionization structure, apart from a sud- 
den chan ge in ionization arou nd spectral type B 1 (the bistability 
jump, see lLamers et al.| [T995). The ionization thus only plays a 
minor role in the wind dynamics and is usually expr essed by a 
small correction to the a parameter (see lAbbottlll982f) . 

As outlined in Sect. 12. 11 the a parameter is determined nu- 
merically within our models and enters the hydrodynamic solu- 
tion scheme as a depth-dependent quantity. Compared with the 
typical values for OB stars {a ~ 0.6) our present models show re- 
markably low values (a « 0.2). This result is in line with our pre- 
vious finding of a » for th e extremely thick winds of WC stars 
dGrafener & Hamannll2005L see also the discussion therein). It 
seems that a — > for thick winds. Such a behaviour is expected 
for the deep atmospheric layers, where the diffusion approxima- 
tion is applicable. For the outer wind such low values again in- 
dicate important qualitative differences to 'classical' line-driven 
winds. At the present stage we attribute this effect to two mech- 
anisms, which are expected for the case of WR-type winds, and 
which are both detectable in our models: the dominance of weak 
spectral lines, and massive line overlap. 
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First, the high density of WR winds leads to an over- 
population of excited energy levels by recombination cascades. 
This results in an increase of the number of optically thin 
lines. Note that a also reflects the relative contribution of op- 
tically thick spe ctral lines to the overall wind driving (see, e.g. 
iPuls et aill2000l) . A dominance of optically thin lines thus leads 
to a decrease of a, or in other words to the formation of a pseudo- 
continuum of thin lines. Second, the resulting ionization struc- 
ture leads to non-local line overlaps. In such a situation it can 
happen that an additional Doppler-shift causes a reduction of the 
radiative force because of the increased line-shadowing. Both 
effects result in a dominance of ionization effects over the usual 
CAK wind physics, as is observed in our models. 

Alternative theories fo r optically thick st ellar winds 
have been proposed by Pi stinner & Eichlerl ( I 19951) and 
iNugis & Lamersl d2002l) . Both works assume that the conditions 
at the critical point depend on the Rosseland mean opacity, and 
they show that the observed mass loss rates of WR stars are in 
agreement with this assumption. In the following we will show 
that the scaling relations obtained from our present models can 
be understood within this framework. They reflect the dominant 
role of the wind optical depth for WR-type mass loss, which we 
identify as the basic reason for the qualitative differences to OB 
star winds. Note, however, that the conditions within our mod- 
els do not strictly reflect the limit of optically thick winds. Our 
present WNL models show only moderate optical depths, so that 
their mass loss actually reflects a case in between the optically 
thick and the optically thin limit. 

The basic assumption of the optically thick wind approach is 
that the radiation field at the critical po int of the wind flow ca n 
be described in the diffusion limit (see lNugis & Lamersll2002l) . 
It follows that 1) the sonic point becomes the critical point of the 
equation of motion u c = a(T c ), 2) a specific (critical) value of 
the Rosseland mean opacity must be reached at the sonic point 
Xc = 4-ncGM+fLi,, and 3) the mass loss rate is connected to 
the density and the temperature at this point via the equation of 
continuity M = 4nR+p c v c . 

The importance of the L/M ratio for this kind of mass loss 
follows directly from condition 2), which means that the radia- 
tive acceleration just balances the gravitational attraction at the 
critical point. For the formation of an optically thick wind it is 
thus mandatory that the L/M ratio and^ c are large enough, i.e., 
the star must be close to the Eddington limit and the mean opac- 
ity at the critical point must be large. Furthermore, the mean 
opacity must increase outward to obtain an outward-directed net 
acceleration above the critical point. These conditions restrict the 
temperature T c at the critical point to limited temperature ranges 
connected to the occurrence of the two Fe opacity-peaks. With 
T c = 30-45 kK our present WNL models are driven by the cool 
Fe-peak opacities. Our previous WCE model with T c = 199 kK, 
on the ot her hand, clearly belonged t o the regime of the hot Fe- 
peak (see Grafener & Hamannll2005l) . 

The necessity to reach certain values of T c marks a basic 
change of paradigm with respect to OB star winds. The temper- 
ature at large optical depth scales as T 4 oc T\t. Hence large wind 
optical depths r c are needed to reach large values of T c . For our 
case, where a specific temperature T c must be reached at r c this 
means that t c oc T~ 4 . 

We found a very strong dependence on for our models 
with fixed T e in Sect. 13. II This behaviour is related to the scal- 
ing properties of t c . If we assume that r c scales with p c and 
the density scale height H, we get r c oc p c H oc p c R+/M e g oc 
p c R+/M. From this it follows directly that t c oc M/M because 
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Fig. 9. WR22: wind acceleration in units of the local gravity 
g. The total atomic density is given as depth index. The blue 
curves indicate the acceleration within our WR22 model from 
Sect. 13. 21 Indicated are the radiative acceleration on free elec- 
trons (a-rhX on continua (a C ontX the total radiative acceleration 
including lines (a rac j), and the total wind acceleration including 
the pressure term (a w ind)- The red line indicates the radiative ac- 
celeration within our test model with a reduced Doppler broad- 
ening velocity vu = 50km/s. 

M = 4nR\p c v c oc R\p c , for the given value of u c = a(T c ). Finally 
we get T+ 4 oc t c oc M/M, or M oc T~ 4 L for fixed L/M. 

This scaling relation indeed predicts a strong dependence on 
T+ and a relatively weak dependence on L for optically thick 
winds. It is caused by the scaling properties of t c with the den- 
sity scale height H, and the requirement that a certain temper- 
ature regime is reached at the critical point. The formation of 
optically thick winds is thus strongly related to the extension of 
the deep atmospheric layers (i.e., the increase of H) close to the 
Eddington limit. 

With M oc T+ 3 5 L 0A our computational models show a simi- 
lar dependence on and an even weaker dependence on L. The 
differences are presumably due to the idealized assumptions in 
our analytic derivation. Nevertheless, we conclude that our mod- 
els show a qualitative behaviour that is very different from OB 
star winds, but can be understood in the framework of the op- 
tically thick wind approach. Note that the specific properties of 
optically thick stellar winds originate from the deep atmospheric 
layers, and that an exact modelling of these layers is required to 
obtain qualitatively correct results. 

4.2. The influence of wind clumping and micro-turbulence 

Apart from the work by ISchmutzl (1 19971) . the influence of wind 
clumping has been largely ignored in previous studies of radia- 
tively driven winds. In our own work we have demonstrated 
that c lumping increases the radiative acceleration in O super- 
giant dGrafener et alj|2002t iGrafener & Hamannll2003h and WR 
star winds ( Grafener & Hamann 120051) . The specific choice of 
the clumping factor D(r) thus has an influence on our results. In 
Sect. |2.2| we have described how D(r) is parametrized within our 
models. This parametrization is based on observational grounds. 
Clumping diagnostics are however very rough and restricted to 
the formation region of the diagnostic features. In principle, D(r) 
thus could be varied in a certain range. In the present section we 
discuss to what extent this would affect our results. In addition, 
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we investigate the influence of the Doppler broadening velocity 
vu- In our models uo is usually set to a fixed, highly supersonic 
value of ud = 100 km/s. This choice is motivated by the observed 
broad ening velocities in P-Cygni type line profiles (e.g. lHamannl 
1 198 lb . and by the fact that our CMF radiative transfer is limited 
to relatively large values of l>d (see Sect. l2.2b . Again, we have a 
certain freedom in the choice of and we want to quantify its 
influence on our results. 

We start with Fig. [9] where we show the wind acceleration 
within our model for WR 22. In this logarithmic diagram we plot 
the wind acceleration a w ; n d in units of the local gravity g, vs. the 
total atomic number density « tot as depth index. The plot covers 
the whole model atmosphere from the hydrostatic layers where 
log(fl w i n d/g) = 0, through the critical point, to the outer boundary 
at r = 1000 R & . The critical point (with v c = 23.3 km/s) and the 
sonic point (with u s = 20.8 km/s) are both located very close to- 
gether at the point where a ra( j ~ g, in agreement with our condi- 
tion 2) from the previous section. In the hydrostatic layers below 
the sonic point (at densities above log(« tot /cirr 3 ) » 12.5), the ra- 
diative acceleration a ra( j (indicated by the dashed blue line) lies 
only slightly below g, with log(a ra d/g) ~ -0.05. This shows that, 
although with F e = 0.55 we are still away from the Eddington 
limit, our model atmosphere is in fact rather close to instabil- 
ity (with a 'true' Eddington parameter, including all opacity 
sources, of T * 0.9). 

Being a late spectral subtype, WR22 is in a temperature 
regime where the cool Fe-peak opacities dominate. In Fig. [9] 
however, we can see that close to the critical point a; 70 % of 
the radiative force is due to continua. Nevertheless, the Fe line 
opacities play the dominant role in determining the mass loss be- 
cause they are responsible for the increase of a ra d which finally 
leads to the crossing of the Eddington limit at the critical point. 
The fact that the line contribution in the deep layers is so small 
compared to continua (which are dominated by Thomson scat- 
tering) indicates that mostly weak lines are responsible for the 
radiative driving, i.e., the sub-critical layers are supported by a 
pseudo-continuum of millions of weak Fe-group transitions. 

From Fig. [9] it also becomes clear why our models are ex- 
tremely sensitive to changes of a ra d close to the critical point. 
The dashed red line indicates the radiative acceleration within a 
model with exactly the same atmospheric structure but with a re- 
duced broadening velocity of i>d = 50 km/s. This causes a rather 
small reduction of a ra d (by roughly -0.025 dex) in the deep at- 
mospheric layers. In the outer wind region the relative effect is 
somewhat larger (-0. 1 dex), the absolute wind velocity is how- 
ever hardly affected because the wind is already close to its ter- 
minal speed. Although the effects in the deep layers are small 
it is obvious that even this small reduction of a la d might cause 
a significant shift of the critical point towards smaller densities 
(i.e., smaller M). A corresponding adjustment of M*, however, 
brings the star back to the same distance from the Eddington 
limit, and should thus restore the previous value of M. 

Test computations indeed confirm this behaviour. In Table[3] 
we compare our WR22 model with full hydrodynamic mod- 
els for which uo has been reduced. For model l>di all parame- 
ters apart from l>d are kept fixed, and M is indeed reduced by 
0.24 dex. For model voi we have compensated this effect by a 
reduction of the stellar mass (0.026 dex) which indeed brings 
the star back to the previous state. Our results do depend on the 
assumed microphysics, because the location of the critical point 
strongly depends on the ratio a ra d/g in the deep atmospheric lay- 
ers. For the same reason our predicted mass loss rates depend 
strongly on the Eddington factor r e . The relatively large value of 
Do = 100 km/s that we have chosen for our computations might 



11 



Model 


WR22 




fD2 


Dl 


D2 




78.1 


78.1 


73.5 


78.1 


78.1 


r c 


0.55 


0.55 


0.585 


0.55 


0.55 


v D / km s -1 


100 


50 


50 


100 


100 


^max 


10 


10 


10 


40 


40 


Tl 


0. / 


0. / 


0. / 


0. / 


U. / 


T2 


0.35 


0.35 


0.35 


0.35 


0.175 


Ml M yr' 


10 -4.85 


10 -4.99 


10 -4.87 


10 -4.61 


10 -4.80 


!)„/ km s" 1 


974 


932 


938 


1859 


1764 




21.7 


25.9 


21.7 


14.5 


18.6 


n 


0.34 


0.23 


0.31 


1.12 


0.69 



Table 3. Model parameters and resulting wind parameters for 
our test models from Sect. 14.21 compared to the WR 22 model 
from Sect. 13. 21 For models i>di and vm the Doppler broadening 
velocity has been reduced to l>d =50 km/s. For models Dl and 
D2 the wind clumping has been increased by a factor of four. 



thus have a significant effect on M. However, it only causes a 
systematic shift of our results towards lower stellar masses, by 
~ 5%. Because the actual conditions in the sub-photospheric 
layers of WR stars are unknown, it is difficult to say if our as- 
sumption is correct. 

The influence of wind clumping turns out to be even more 
important. In the last two columns of Table|3]we summarize the 
results of two test computations where the maximum clump- 
ing factor D max has been scaled up by a factor of 4. Note that 
the spectroscopic properties of WR stars are determined by the 
product M \[D (for details, see lHamann & Koesterkell99~8"ah . An 
increase of the clumping factor by a factor of four thus approxi- 
mately mimics the opacity of a wind with two times higher den- 
sity. Consequently, a la d also should scale with Vd. Note that 
such a scaling relation is only expected to hold strictly for the 
populations of excited energy levels which are dominated by re- 
combination processes. 

Our test computations (models Dl and D2 in Table|3) show 
that the expected scaling of a ra d °c Vd is fulfilled with an aston- 
ishing precision. Both test models with £> max = 40 show nearly 
exactly twice the terminal wind speed as our Z) max = 10 mod- 
els. Again, this is a hint that weak line transitions from excited 
(recombination-dominated) levels play a major role in the wind 
driving of WR stars. Moreover, the radial dependence of the 
wind clumping is also important. With our standard prescription 
of a r-dependent clumping factor D, which rises from D = 1 to 
Anax between Rosseland optical depth t\ = 0.7 and T2 = 0.35 
(see Sect. |2.2| >, the increase of D max also leads to a significant 
increase of M (model Dl in Table|3]l. The reason for this is that 
in our prescription of D(t) the wind clumping already sets in be- 
low the critical point, and thus affects the critical condition. For 
model D2 we have set T2 to a value of 0. 175. In this way we have 
modified the growth rate of D(r) such that the conditions at the 
critical point are similar to our WR22 model. Again, the mass 
loss reacts very sensitively to the change of D(r) and returns to 
the same value as the WR 22 model. 

We conclude that different assumptions in the microphysics 
of our models may indeed cause significant changes in the pre- 
dicted mass loss rates, if the region around the critical point is 
affected. However, these effects can always be compensated by a 
slight re-adjustment of the prescribed stellar mass. For the wind 
clumping we find that a ra d scales with Vd. An increase of D in 
the outer wind layers only leads to an increase of the terminal 
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Fig. 10. The Z-dependence of WNL star mass loss: comparison 
of different mass loss predictions. The solid blue lines indicate 
the Z-dependent mass loss rates according to Eqs. (ID)-© f° r 
fixed values of F e = 0.3-0.9 (note that only the thick parts of the 
blue curves are in the region where grid models have actually 
be en computed). The y are compared to the mass loss predictions 
by IVink et al.l (1200 ll solid red lines) for T e = 0.3, 0.6, and 0.9 
and th e mass loss relation for WNL stars by IVink & de Koterl 
(I2005L dashed red line), which was obtained with F e = 0.34. 



wind speed by a factor Vd. The possible discrepancies of our 
predicted wind speeds with observations (Sect. l3.3t thus might 
be resolved by a strong increase of D(r) in the outer wind. 

4.3. Comparison with other mass loss predictions 

In the last few years an impressive agreement between theoreti- 
cal mass loss predictions for OB stars and observations has been 
achieved (but s ee the discussion in Sec t. 5.31>. On the theoretica l 
side models bvlVink et al.ld2000l l200ll). Pauldrach et al. (2001). 
and iKrticka & Kubatl d20O4l) agree well. Observationally, these 
results are confirmed on a large scale by the VLT-F LAMES sur- 
vey of massive stars dMokiem et al.l 120051 l2007al) . For the Z- 
dependence an exponential relatio n of the form M oc Z 7 is gen- 
erally accepted. I Vink et all d2001l) find an exponent of y « 0.85 
for a broad range of OB star parameters. Note that this exponent 
is only valid if v m is kept fixed. It thus needs to be corrected 
for the dep endence of Voo(Z ), resulting in a relation of the form 
M oc Z 69 dVinket al.ll200ll) . 

For WR stars the situation is less clear be cause only a few 
model s are available. For late spectral subtvpes lVink & de Koterl 
d2005l) recently claimed a similar Z-dependence as for OB stars, 
with y = 0.86. As discussed previously, our present WNL mod- 
els show a qualitatively different behaviour. Apart from the in- 
fluence of clumping, micro-turbulence, and high wind density, 
it seems that models using a detailed CMF radiative transfer 
give results different from wind models based on the Sobolev ap- 
proximation. Without wind clumping both available CMF codes 
(CMFGEN bv lHillier & Miiieil (U9981) and our own PoWR code) 
tend to predict very low wind accelerations a ra d- The models 
can only account for the observed wind momenta if clump- 
ing factors of the order of D = 10 or larger are assumed 



dGrafeneret alfcOOOtlGrafener et alJl2002tlGrafener & Hamannl 
120051: iHerald et alj|2001h lHilUeill2003h iHillier et al.ll2003l) . The 
re ason for this is not yet fully understood (but see the discussion 
in iHillier'e t al. 2003]). In the following we try to quantify these 
differences by a detailed comparison with previous works. 



For comparison we use the Z - depen dent mass loss prescrip- 
tio n for OB stars by IVink et al.l d2001l) . and the WNL models 
by I Vink & de Koterl d2005l) . Note that, although the modeli ng of 
WNL subtypes was not particularly intended bv lVink et al.l their 
model grid partly overlaps with our present computations. In 
Fig.[l0]we have extracted the Z-dependence of our WNL models 
from Eqs. (H)-©, for fixed values of T e . The blue curves indicate 
our results for Eddington parameters F e =0.3 to 0.9 (note that 
only the bold parts of the curves belong to the parameter range 
where models have actually been computed). The stellar param- 
eters are fixed to log(r+/K) = 4.65, log(L+/L ) = 6.3, and 
Xh = 0.4. The r ed solid lines indi cate the mass loss predictions 
for OB stars by IVink et al.l (1200 ll) for the same stellar parame- 
ters and the same range of Eddingto n factors. The dashed re d 
curve represents the WNL models by I Vink & de Koterl d2005l) . 
for slightly different stellar parameters (see b elow). 

The OB star models bv lVink et al.l (1200 ll) have a strictly ex- 
ponential Z-dependence. Although F e spans the same parameter 
range as our WNL models, the resulting mass loss rates cover a 
much narrower range, i.e., the dependence on F e is much weaker. 
For large Eddington factors we obtain much higher mass loss 
rates, and a flatter decline with Z. These important differences 
are most likely cause d by the diff erent treatment of the wind dy- 
namics in both codes. I Vink et al.l prescribe the velocity structure 
in their models and solve for a global consistency of the wind en- 
ergy. The peculiar wind physics in the deep layers of WR winds, 
as described in Sect. 14. II is not considered in such an approach. 

The WNL models by I Vink & de Koterl d2005l) are indicated 
by the red dashed line in Fig. [10] They are computed for a sim- 
ilar effective temperature as our models (log(T + /K) = 4.60), 
but for a lower hydrogen surface abundance (Xu = 0.15), and 
a much lower luminosity (log(L*/L ) = 5.62). The stellar mass 
is fixed to 20 M , corresponding to F e = 0.37, and the terminal 
wind speed is set to 840 km s -1 . For Z < Z these models show 
exactly the same Z-dependence as their OB star counterparts, 
while for Z > Z the mass loss tends to saturate. Moreover, 
IVink & de Koterl detect a flattening of the Z-dependence for 
Z < 1O~ 4 Z , which is beyond the parameter range investigated 
here. 

For Z > Z > 10~ 4 Z Q the WNL models bv lVink & de~Koterl 
(120051) are i n exact agre e ment w ith the ma ss loss prescription for 
OB stars bv lVinketal] d200ll) . Note that I Vink et all d200ll) as- 
sume a ratio of u<x,/u e sc = 2.6 for OB stars. With a typical stellar 
mass of 43.5 M for an OB star of this luminosity their mass loss 
prescription gives M = 10~ 569 M yr~' and = 2520 km s~' at 
solar metallicity. A reduction of Voo to 840kms~ I leads to an in- 
crease of the mass loss rate by a factor of four, to 10~ 510 M Q yr~' . 
A reduction of the stellar mass to 20 M together with an ad- 
justment of the Eddington factor to T e = 0.37 further increases 
M by a factor of 1.5, to M = 10' 491 M yr'. This is in ex- 
act agreement with the WNL model by I Vink & de Koterl (120051) 
with M = 10~ 489 M yr -1 . Because the Z-dependence also is the 
same, the model s bvlVink et al.l (1200 1[) th us closely resemble the 
WNL models bv lVink & de Koterl d2005l) . 

According to Eqs. (|3]l-(|6]l our own models give a much lower 
value of M = 10~ 607 M Q yr~' . Note, however, that our formula is 
not reliably applicable because the desired Eddington factor of 
0.37 is very close to the parameter Fo(Z ) = 0.33 (see the discus- 
sion at the end of Sect. 13. A\ . A direct model computation yields 
M = lO _5 4 M yr _I (see Table|2]for details). The discrepancy is 
caused by the influence of the He n continuum, which effectively 
increases the Eddington factor close to the critical point in this 
specific parameter range. Neverth eless, our mass los s rate is still 
significantly below the value from I Vink & de Koterl 
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We conclude that our models are much more sensitive to the 
adopted Eddington factor than others. This leads to higher or 
lower mass loss rates, dependent on F e . According to our mod- 
els the proximity to the Eddington limit is the primary reason 
for the enhanced WR-type mass los s . This possibility has also 
been discussed by IVink & de Kot er (2005). Our models how- 
ever show qualitative changes in the optically thick regime that 
are also reflected in the Z-dependence of the mass loss. While 
the other studies find a universal relation with M oc Z° 85 for al- 
most all types of hot stars, we detect a flattening of the M(Z) 
relation for F e > 0.7. For lower values of T e our models show 
a similar Z-dependence to that determined by IVink & de Koterl 
The important fact that our models yield high mass loss rates 
over the whole range of metallicities implies that WR-type mass 
loss might gain importance in metal-poor environments (see the 
discussion in Sect.l - 

5. Discussion 

In the preceding sections we have presented the first fully self- 
consistent models for the winds of WNL stars. These models 
reveal important information about the nature of these objects 
and their mass loss. In the following we discuss the most impor- 
tant points: the reason for their enhanced mass loss (Sect. 15. it , 
their evolutionary state (Sect. l5.2t , the effect of wind clumping 
(Sect. l5.3l l, the role of radiative wind driving (Sect. l5.4l ). the Z- 
dependence of WR-type mass loss (Sect. 15.51 ), and its potential 
role in extremely metal-poor environments (Sect. |5.6l l. 

5.1. What is the reason for WR-type mass loss? 

The most important conclusion from our present work is that 
WR-type mass loss is primarily triggered by the proximity to 
the Eddington limit. As the main reason we identify the exten- 
sion of the deep atmospheric layers, leading to the formation of 
optically thick stellar winds. We have performed the first detailed 
computations of the atmospheric structure of such objects, and 
we have identified qualitative differences to OB star winds. The 
most important is the strong dependence of the mass loss rates 
on the Eddington factor F e and on the effective stellar temper- 
ature r*. The idea that stellar mass loss increases close to the 
Eddingto n limit is of cour s e not new. For example, for the case 
of LBVs IVink & de Koterl d2002l) also detected a strong depen- 
dence on F e . Nevertheless, in stellar evolution models the surface 
composition and temperature are commonly used to discriminate 
between different states of mass loss. 

The proximity to the Eddington limit provides a natural ex- 
planation for the occurrence of the WR phenomenon, for He- 
burning stars with their inherently high L/M ratios as well as for 
over-luminous H-burning stars. The latter is expected for very 
massive stars at the end of their main-sequence evolution and 
for f ast rotators where fresh h ydrogen is mixed into the stellar 
core dMevnet & Maederll2000l) . Note that such objects are also 
expected to show WN surface compositions due to rotationally 
induced mixing and mass loss. 

The wind driving in our present models mainly originates 
from line opacities of Fe m-vm. These ions belong to the "cool" 
Fe opacity peak, in contrast to our previous WCE models, where 
the "hot" F e peak (Fe ix-xvi) is respon sible for the driving of the 
wind base ( Grafener & Hamannl2005l) . We thus confirm the pre- 
diction bv lNugis & Lamersr ( l2002l) that the winds of WNL stars 
are driven by the cool Fe peak opacities, whereas for early spec- 
tral subtypes the critical point is located in much deeper layers 
where the hot Fe peak dominates. 



5.2. The evolutionary status of WNL stars 

The dichotomy among the WN subtypes, where the H-rich 
WNL stars show much higher lu minosities than the earlier sub- 
types (e.g., lHamann et al.ll2.006h . indicates that the progenitors 
of H-rich WNL stars are very massive stars with masses around 
120 M & . In Sect. 13. II we have shown for the example of WR22 
(WN 7h) that the observed spectrum, including the absolute flux 
distribution, is reproduced by a very luminous model (L+ = 
10 63 L o ) with T e = 0.55, corresponding to a stellar mass of 
78 M e . Such a high mass implies that WR 22 is still in the phase 
of central H-burning. Our models thus suggest that, for very 
massive stars, the WNL phase might already occur at the end 
of the main-sequence evolution, when the L/M ratio increases 
due to the increase of the mean molecular weight in the stellar 



Our parameters for WR 22 are in agreem ent with the evolu- 
tionary model for a rotating 120 M Q star by iMevnet & Maederl 
(120031) . which is still in the phase of central H-burning (at an 
age of 2.2 10 6 yr). This means that our derived mass suggests 
that WR22 is in the end phase of the main-sequence evolu- 
tion. The observed mass loss rate of WR22 (10~ 4 85 M G yr~') 
is not in agreement with the evolutionary models. In the mod- 
els the mass loss switches from O-star mass loss (in the range 
of lO^-lO^Moyr 1 ) to WNL mass loss (~ lO^M^r 1 ) 
when the hydrogen surface mass fraction falls below 0.4. This 
sudden switch reduces the stellar mass very rapidly. The pre- 
vious evolution of WR 22 has probably been different because, 
even with its present mass loss rate as the average value, it can 
only have lost up to 28 M Q within 2 Myr. 

Our models suggest an H-rich WNL phase at the end of the 
main sequence evolution of very massive stars. Within the stan- 
dard evolutionary scenario such o bjects enter the L BV stage af- 
ter leaving the main sequence (e.g.. ICrowtherl2007l) . After losing 
their H-rich envelopes, they presumably evolve into the 'classi- 
cal' He-burning WR stages. Our models thus support an evolu- 
tionary sequence of the form O — > WNL (H-rich) — > LBV — > 
WN(H-poor) — » W C for very massive stars, as proposed by 
lLanger et al.l ([1994). According to our results in Sect. 13. II the 
mass loss from WNL stars increases strongly for decreasing stel- 
lar temperatures (M oc r + 3 5 ). For temperatures below the values 
examined here, this might indeed lead to a smooth transition into 
the quiet LBV phase. 

Additional evidence that the most massive stars are WNL 
stars comes from observations of young stellar clusters. There 
are several cases where the brightest stars in such clusters are 
WNL stars. A good example is the Carina OB cluster Trl6 
where only two evolved stars are prese nt, the extremely lu- 
minous WNL star WR25 (WN6h, see lHamann et ail 120061) 
and 77 Car. Also t he br i ght WNL stars in t he Arches clus- 
ter jNajOTojetallj2003 IMEnnsitiD llQQl) and in R136a 



(Ide Koter et al.l 19971) fall into the same category as WR22 
and WR 25 - extremely luminous H-rich WNL stars with rel- 
atively weak emission lines. Moreover, the hi ghest masses in bi- 
nary systems are measured for WNL stars (iRauw et ail 19961: 
I Schweickhardt et al.l 1 19991: iBonanos et ail 120041: iMoffat et all 
|2007|) . H-rich WNL stars thus probably mark an important phase 
of strong mass loss at the end of the main-sequence evolution of 
very massive stars. 

Note that we do not exclude the possibility of less luminous 
WNL stars being in the phase of central He-burning, succeeding 
a RSG/LBV phase. Our present assumption of high luminosities 
relies on galactic WNL stars with known distance, which are 
usually members of dense clusters. There are, however, many 
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WNL stars that are not located in clusters, and that might have 
much lower luminosities. Such a case is WR40 (WN 8h) which 
is located in the center of the ring nebula R CW 58, a possi ble 
remnant of a previous RSG outflow (see, e.g.. IMarstonll 1 9991) . 

5.3. The effect of wind clumping 

Wind clumping is routinely taken into account in spectral analy- 
ses of WR stars. Its main effect is that the mean (p 2 ) increases 
and thus the derived mass loss rates, relying on p 2 -diagnostics, 
decrease with M oc 1/ Vd (see lHamann & Koesterkdll998ah . 
In the present work we have shown that clumping also affects 
the wind acceleration. The reason for this is that excited energy 
levels are more populated by recombination processes. For the 
dense winds of WR stars, the resulting increase of optically thin 
lines leads to an increase of the radiative acceleration, approxi- 
mately with a ra d oc Vd. Moreover, the force multiplier param- 
eter a is reduced. At this point it is necessary to emphasize that 
we assume small-scale clumps (micro-clumping) in our models. 
In this limit the separation between clumps is so small that the 
mean opacity can be used in the radiative transfer. In the limit 
of large clump separation (macro-clu mping) the opt i cal de pth of 
individual clumps may become large. IBrown et al.l ((2004) have 
shown that in this case the radiative acceleration is reduced by 
ge ometrical effects . With a modified approach for the wind opac- 
ity lOskinova et al.l d2007l) have shown that spectral diagnostics 
also are considerably affected by macro-clumping. 

We thus have two counteracting effects that may affect our 
results. Micro-clumping increases the radiative acceleration with 
flrad K VD, and macro-clumping decreases a ra d dependent on 
the detailed clump-geometry. In our present models only micro- 
clumping is taken into account. We have shown that it only af- 
fects the wind speed as long as it is restricted to the outer wind 
region above the critical point. In some of our models the criti- 
cal region is however marginally affected by clumping. This of 
course introduces an uncertainty in our mass loss rates. 

Note that significant clumping factors in the range of 
D — 10-100 have recently been proposed for OB star winds 
(ICrowther et all l200l iHillier et alJ l200l iBouret et al.l l200l 
2005; Ful lerton et al T l2006h . If the effects of micro-clumping 
dominate, the observed mass loss rates of OB stars thus should 
be revised downward by factors of 3-10. This would of course 
introduce a discrepancy in the previous wind models for OB 
stars because these do not take clumping into account. As we 
discussed in Sect. 14. 31 our models yield a much lower radia- 
tive force than others, which seems to be connected to the de- 
tailed CMF radiative transfer. In contrast to the standard Sobolev 
models, spec tral lines ha ve a finite width in our approach. 
Interestingly, iLucvl d2007l) has proposed that the assumption of 
infinitely narrow spectral lines leads to an over-estimation of 
the mass loss rates in the standard models, i.e., OB star mass 
loss rates should actually be lower. Because of its influence on 
the wind dynamics and the wind diagnostics, it thus seems that 
clumping is fundamentally important for the physics of radia- 
tively driven winds. Because the actual mechanisms of wind 
clumping are presently not well understood, it introduces sig- 
nificant uncertainties in all types of mass loss predictions for hot 
stars. 



5.4. Are WR winds generally driven by radiation? 

In the present work we have shown that WNL star winds can 
be explained self-consistently by pure radiative driving. As the 



reason for their strong optically thick winds we have identi- 
fied the extension of the deep atmospheric layers due to the 
radiative force on Fe line opacities. For our present WNL 
models the cool Fe-peak opacities are responsible for this ef- 
fect. Together with our prev ious results for early-type WC stars 
JGrafener & HamannI 120051) . where the inner wind is driven 
by the hot Fe-peak, one might now argue that WR winds are 
generally driven by radiation. This would imply a general Z- 
dependence of WR mass loss. 

The fact that the two Fe opacity peaks act in different tem- 
perature regimes, however, has important consequences for our 
models. First, because of the gap between the two peaks, the 
mass loss cannot be maintained for intermediate stellar temper- 
atures (T* = 60-100 kK). Second, for cool objects, the hot Fe- 
peak disturbs the inner atmospheric structure. 

The first point is connected to the well-known "extension 
problem" of H-free WR stars. For H-free WN and WC stars 
stellar structure models predict much smaller radii than ob- 
served. Note that the predicted radii correspond to stellar tem- 
peratures above ~ 100 kK, in agreement with our wind models. 
Obs erved stellar tempera tures, however, extend down to ~ 50 kK 
(e.g. lHamann et all 2006). From a theoretical point of view such 
objects are thus not expected to exist, and they are not expected 
to have WR-type stellar winds. Ilshii et all d 19991) have shown 
that very luminous WR stars might form an extended convection 
zone close to the stellar surfac e which can account for the ob- 
served extension. According to lPetrovic et al.l d2006l) this exten- 
sion however only occurs under the assumption of hydrostatic 
equilibrium. As soon as dynamic terms are included in their 
computations the extension disappears. Interestingly, this exten- 
sion is caused by the hot Fe-peak. According to our computa- 
tions the hot Fe-peak drives the stellar wind, i.e., the possible 
extension of the outer layers of WR stars is hindered by their 
mass loss. A possible solution for this problem could be the pul- 
sational driving of the deep layers (e.g., by strange mode insta- 
bilities as proposed bv lGlatzel & Kaltschmidtll2002l) . This could 
help to overcome the opacity-gap between the two Fe-peaks, 
and also account for the observed extension. Another possibility 
would be to bring the stars closer to the Eddington limit. For He- 
bu rning stars the L/M ratio is however fixed. As also discussed 
bv lVink & de Koteri d2005l) . fast stellar rotation may help to re- 
duce the surface gravity of such objects and initiate strong mass 
loss. That this process works for optically thick winds, however , 
remains to be shown. For OB star winds Madura et al. d2007l) 
have demonstrated that rotation increases the mass loss at most 
by a factor of two. 

The second problem only arises for late spectral subtypes. 
For these stars the cool Fe-peak drives the stellar wind, and the 
hot Fe-peak is located in deep layers at temperatures around 
160kK. For a star close to the Eddington limit it can thus hap- 
pen that the radiative force exceeds gravity in the deep hydro- 
static layers. In the stellar interior this would lead to the onset of 
convection. Particularly for stars with dense winds this situation 
already occurs at the inner boundary of our atmosphere mod- 
els. We thus reach a mass loss limit above which we cannot ob- 
tain stationary wind mo dels. In reali t y this presumably leads to a 
non-stationary situation. iDorfi et all d2006l) have recently shown 
that the observed variability in WN 8 subtypes can indeed be ex- 
plained by instabilities that are caused by the hot Fe-peak. 

The comparison with observations in Sect. 13.11 has shown 
that only a part of the galactic WR sample is reproduced by our 
models. The WN 8 subtypes particularly have higher wind densi- 
ties than predicted. The comparison with the WN 8 model from 
IVink & de Kotej d2005l) in Sect.gJ] has shown that this prob- 
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lem still holds if lower luminosities are adopted. According to 
our models a different driving mechanism thus would be de- 
sirable, at least for the wind base of these objects. Again, this 
mechanism could be related to the observe d pulsations (see e.g. 
iMarchenko et alj!998l:lL~efevre et al.l2005h or rotational driving. 
Note that also the non-detection of X-ray emission from WN 8 
subt ypes might be a hint to a different wind driving mechanism 
(see IOskinovall2005l) . 



5.5. The Z-dependence of WR mass loss 

As expected for radiatively driven winds, our models show a 
strong Z-dependence. Because of the important influence of the 
Eddington factor, it is however not sufficient to parametrize the 
mass loss as a function of Z alone. In Sect. 13. 41 we have thus 
prepared a mass loss recipe which takes the combined effect of 
T e and Z into account. According to Eqs. (IUi-© the metallic - 
ity Z affects the resulting mass loss rates in two ways. First, the 
parameter r (Z) decreases for increasing Z (Eq.[6]). To character- 
izes the limiting Eddington factor for which WR-type mass loss 
sets in. An increasing metallicity thus brings the star effectively 
closer to the Eddington limit, and supports the formation of WR- 
type winds. Second, the exponent /3(Z) increases with Z (Eq.[5]). 
This results in the flattening of the M(Z) relation for low Z and 
large T e (see Sect. |4.3t . 

A higher metallicity thus has two effects. The mass loss 
from WR stars increases and WR stars are formed more eas- 
ily. According to our models, at solar metallicity the most mas- 
sive stars already enter the WNL phase at the end of the In- 
huming stage. For Z = 3Z we could even show that stars 
with 120 M H already become WNL s tars on the zero-age main 
sequence dGrafener & Hamamill2006l) . In low metallicity envi- 
ronments, on the other hand, the occurrence of WR-type winds 
should be restricted to He-burning objects and maybe to fast ro- 
tators. The influence of metallicity on WR populations is thus 
extremely complex. In different environments different objects 
become WR stars. 



5.6. WR-type mass loss in extremely metal-poor 
environments 



In Sect. l4.3l we have shown that the slope of the M(Z) relation 
becomes shallower for increasing F e . Close to the Eddington 
limit our models are thus able to produce large mass loss rates 
also at very low Z. Note that these are presently the only wind 
models for low metallicities for which such an efficient mass 
loss mechanism has been detected. The large Eddington factors 
that are required (F e ~ 0.8) are however only reached by very 
massive objects or stars close to critical rotation. Both types of 
o bjects are indeed expec ted at low Z (see also the discussion 
in lvink & de Kotei1l2005l) . For extremely low met allicities (Z < 
10~ 4 Zp) more massive stars are expec ted to form dBromm et alj 
119991: iNakamura & Umemural 120021) . Moreover, evolutionary 
calculations indicate that rotating mas sive stars tend to approac h 
critical rotation more easily at low Z dMevnet & Mae der 2002). 
Our present models show that high mass loss rates can be main- 
tained in this regime, and that such objects should appear as 
WR stars with particularly low wind velocities. Note, however, 
that the question of whether rotation really leads to increased 
mass loss r emains to be investig ated in more detail. For CAK- 
type winds iMadura et al.l d2007 ) recently demonstrated that ro- 
tation increases M at most by a factor of two. 



A new and very important result from our calculations con- 
cerns the role of primary elements for the wind d riving of ex- 
tremely meta l -poor stars. lYoon & Langerl d2005l) . as well as 
iMevnet et al.l d2006l) . have shown that rotating massive stars at 
extremely low Z (1O~ 5 Z and 1O _8 Z respectively) are capable 
of producing nitrogen mass fractions in the range of one percent 
at their surface. The mass loss from such objects may potentially 
enrich the early ISM with freshly produced ni t rogen. With the 
adopte d stellar yields from lMevnet et alJd2006l) . IChiappini et al.l 
d2006l) could explain the observ ed high nitrogen abundances in 
extremely metal-poor halo stars (Spite et al.l 20051). 

Our iron-free models in Sect. l3.3l now confirm that such self- 
enriched objects can develop strong stellar winds. We find that 
iron-free objects with solar Z (i.e., Xqno = Z Q and Xp e = 0) can 
reach mass loss rates similar to stars with standard composition 
at Zq/50. This means that medium-complex elements like C, N, 
and O are much less efficient for the wind driving than the com- 
plex Fe-group elements. For more realistic mass loss predictions 
it is thus important to take the relative contributions of different 
elements into account. 

For the application of Z-dependent mass loss relations like 
our Eqs. @, (0, and (O, we thus suggest the use of an effec- 
tive metallicity Z e g, where the relative contributions of different 
elements are adequately weighted. In the present case of CNO- 
enriched WNL stars we propose an approximation of the form 



Z e ff 



1 ^CNO X 



Fe 



50 XqnO.Q Xp^Q 



(7) 



6. Conclusions 

We have performed a systematic study of the mass loss from 
late-type WN stars by means of a new generation of non-LTE at- 
mosphere models, including a self-consistent treatment of the 
wind hydrodynamics. This was the first application of these 
complex models on a large scale. We could show that the strong 
winds from WNL stars can be explained self-consistently by ra- 
diative driving, and that the proximity to the Eddington limit 
is the primary trigger of the enhanced WR-type mass loss. 
Moreover we have identified important qualitative differences to 
OB star winds that are caused by the altered wind physics. 

Our models reproduce the observed sequence of late spec- 
tral subtypes from WN6 to WN9 qualitatively. The obtained 
wind densities are in agreement with the galactic WNL sample, 
except for the enigmatic WN 8 subtypes with strong winds. A 
detailed comparison with the spectroscopic binary WR 22 sup- 
ports our key statement that H-rich WNL stars are very massive, 
luminous stars close to the Eddington limit. Our models imply 
that these objects are still in the phase of central H-burning, in 
line with an evolutionary sequence of the form O — » WNL (H- 
rich) — » LBV — » WN(H - poor) — » WC fo r very massive stars 
(e.g., lLanger et all Il994t Icrowfheri 120071) . Due to the strong 
temperature-dependence of the mass loss (M oc 7^ 3 ' 5 ), an evolu- 
tion towards cooler effective temperatures may indeed lead to a 
continuous transition into the quiet LBV phase. 

As expected for radiatively driven winds, our models show 
a strong Z-dependence. The influence of the Eddington factor 
r e however is at least equally important. We find that the M(Z) 
relation for WNL stars becomes flatter for higher F e , i.e., high 
mass loss rates can also be maintained for very low Z. Our 
models thus predict an efficient mass loss mechanism for low 
metallicity stars. For stars in extremely metal-poor environments 
(i.e. the second generation of massive stars) we find that the 
surface-enrichment with primary nitrogen, as described, e.g. by 
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iMevnet et all d2006h . may lead to strong mass loss. We thus con- 
firm that these first WN stars might play a key role in the enrich- 
ment of the early ISM with freshly produced nitrogen. 



References 

Abbott, D. C. 1982, ApJ, 259, 282 

Anderson, L. S. 1989, ApJ, 339, 558 

Auer, L. H. & Mihalas, D. 1970, MNRAS, 149, 65 

Barniske, A., Hamann, W.-R., & Grafener, G. 2006, in ASP Conf. Ser. 353: 
Stellar Evolution at Low Metallicity: Mass Loss, Explosions, Cosmology, ed. 
H. J. G. L. M. Lamers, N. Langer, T. Nugis, & K. Annuk, 243 
Bonanos, A. Z., Stanek, K. Z., Udalski, A., et al. 2004, ApJ Lett., 61 1, L33 
Bouret, J.-C, Lanz, T., & Hillier, D. J. 2005, A&A, 438, 301 
Bouret, J.-C, Lanz, T., Hillier, D. J., et al. 2003, ApJ, 595, 1182 
Bromm, V., Coppi, P. S., & Larson, R. B. 1999, ApJ, 527, L5 
Brown, J. C., Cassinelli, J. P., Li, Q., Kholtygin, A. F., & Ignace, R. 2004, A&A, 
426, 323 

Castor, J. I., Abbott, D. C, & Klein, R. I. 1975, ApJ, 195, 157 
Chiappini, C, Hirschi, R., Meynet, G., et al. 2006, A&A, 449, L27 
Conti, P. S., Garmany, C. D., & Massey, P. 1989, ApJ, 341, 113 
Crowther, P. A. 2000, A&A, 356, 191 
Crowther, P. A. 2007, ARA&A, 45, 177 
Crowther, P. A. & Dessart, L. 1998, MNRAS, 296, 622 
Crowther, P. A., Hillier, D. J., Evans, C. J., et al. 2002, ApJ, 579, 774 
Cunto, W. & Mendoza, C. 1992, Rev. Mex. Astron. Astrofis., 23, 107 
Davidson, K., Smith, N, Gull, T. R., Ishibashi, K., & Hillier, D. J. 2001, AJ, 
121, 1569 

de Koter, A., Heap, S. R., & Hubeny, I. 1997, ApJ, 477, 792 
Dorfi, E. A., Gautschy, A., & Saio, H. 2006, A&A, 453, L35 
Fitzpatrick, E. L. 1999, PASP, 1 1 1, 63 

Fullerton, A. W., Massa, D. L., & Prinja, R. K. 2006, ApJ, 637, 1025 

Glatzel, W. & Kaltschmidt, H. O. 2002, MNRAS, 337, 743 

Gosset, E., Remy, M., Manfroid, J., et al. 1991, Informational Bulletin on 

Variable Stars, 3571, 1 
Grafener, G. & Hamann, W.-R. 2003, in IAU Symp., Vol. 212, A Massive 

Star Odyssey: From Main Sequence to Supernova, ed. K. A. van der Hucht, 

A. Herrero, & E. Cesar (San Francisco: ASP), 190 
Grafener, G. & Hamann, W.-R. 2005, A&A, 432, 633 

Grafener, G. & Hamann, W.-R. 2006, astro-ph/0609672 to appear in "The Metal 

Rich Universe" (Ed. G. Israelian) 
Grafener, G, Hamann, W.-R., & Koesterke, L. 2000, in ASP Conf. Ser., Vol. 204, 

Thermal and Ionization Aspects of Flows from Hot Stars, ed. H. J. G. L. M. 

Lamers & A. Sapar (San Francisco: ASP), 215 
Grafener, G., Koesterke, L., & Hamann, W.-R. 2002, A&A, 387, 244 
Hamann, W.-R. 1981, A&A, 93, 353 
Hamann, W.-R. 1985, A&A, 148, 364 
Hamann, W.-R. 1986, A&A, 160, 347 
Hamann, W.-R. & Grafener, G. 2003, A&A, 410, 993 
Hamann, W.-R. & Grafener, G. 2004, A&A, 427, 697 
Hamann, W.-R., Grafener, G., & Liermann, A. 2006, A&A, 457, 1015 
Hamann, W.-R. & Koesterke, L. 1998a, A&A, 335, 1003 
Hamann, W.-R. & Koesterke, L. 1998b, A&A, 333, 251 

Hamann, W.-R., Leuenhagen, U., Koesterke, L., & Wessolowski, U. 1992, A&A, 
255, 200 

Hamann, W.-R., Pefia, M., Grafener, G., & Ruiz, M. T. 2003, A&A, 409, 969 
Herald, J. E., Hillier, D. J., & Schulte-Ladbeck, R. E. 2001, ApJ, 548, 932 
Hillier, D. J. 2003, in IAU Symp., Vol. 212, A Massive Star Odyssey: From Main 

Sequence to Supernova, ed. K. A. van der Hucht, A. Herrero, & E. Cesar (San 

Francisco: ASP), 70 
Hillier, D. J., Lanz, T., Heap, S. R., et al. 2003, ApJ, 588, 1039 
Hillier, D. J. & Miller, D. L. 1998, ApJ, 496, 407 
Howarth, I. D. & Schmutz, W. 1992, A&A, 261, 503 
Hummer, D. G. & Mihalas, D. 1988, ApJ, 331, 794 
Iglesias, C. A. & Rogers, F. J. 1996, ApJ, 464, 943 
Ishii, M., Ueno, M., & Kato, M. 1999, PASJ, 51, 417 
Koesterke, L. & Hamann, W.-R. 1995, A&A, 299, 503 
Koesterke, L., Hamann, W.-R., & Grafener, G. 2002, A&A, 384, 562 
Koesterke, L., Hamann, W.-R., & Kosmol, P. 1992, A&A, 255, 490 
Krticka, J. & Kubat, J. 2004, A&A, 417, 1003 
Kudritzki, R. P., Puis, J., Lennon, D. J., et al. 1999, A&A, 350, 970 
Kurucz, R. L. 1991, in NATO ASI Series, Vol. 341, Stellar Atmospheres - 

Beyond Classical Models, ed. L. Crivellary, I. Hubeny, & D. G. Hummer, 

441 

Kurucz, R. L. 2002, in AIP Conf. Proc. 636: Atomic and Molecular Data and 
Their Applications, 134 



Kurucz, R. L. & Bell, B. 1995, Kurucz CD-ROM No. 23. Cambridge, Mass.: 

Smithsonian Astrophysical Observatory 
Lamers, H. J. G. L. M. & Nugis, T. 2002, A&A, 395, LI 
Lamers, H. J. G. L. M., Snow, T. P., & Lindholm, D. M. 1995, ApJ, 455, 269 
Langer, N, Hamann, W.-R., Lennon, M., et al. 1994, A&A, 290, 819 
Lefevre, L., Marchenko, S. V, Moffat, A. F. J., et al. 2005, ApJ Lett., 634, L109 
Leuenhagen, U. & Hamann, W.-R. 1994, A&A, 283, 567 
Leuenhagen, U., Hamann, W.-R., & Jeffery, C. S. 1996, A&A, 312, 167 
Lucy, L. B. 1964, in First Harvard-Smithsonian Conference on Stellar 

Atmospheres, ed. E. H. Avrett, O. J. Gingerich, & C. A. Whitney, 

Smithsonian Astrophysical Observatory Special Report No. 167 (Cambridge, 

Mass.), 93 
Lucy, L. B. 2007, A&A, 468, 649 
Lucy, L. B. & Abbott, D. C. 1993, ApJ, 405, 738 
Lundstrom, I. & Stenholm, B. 1984, A&AS, 58, 163 
MacFadyen, A. I. & Woosley, S. E. 1999, ApJ, 524, 262 
Madura, T. I., Owocki, S. P., & Feldmeier, A. 2007, ApJ, 660, 687 
Marchenko, S. V, Moffat, A. F. J., Eversberg, T., et al. 1998, MNRAS, 294, 642 
Marston, A. P. 1999, in IAU Symp. 193: Wolf-Rayet Phenomena in Massive 

Stars and Starburst Galaxies, ed. K. A. van der Hucht, G. Koenigsberger, & 

P. R. J. Eenens, 306 
Martins, F., Hillier, D. J., Paumard, T, et al. 2008, A&A, 478, 219 
Massey, P. & Johnson, J. 1993, AJ, 105, 980 
Meynet, G, Ekstrom, S., & Maeder, A. 2006, A&A, 447, 623 
Meynet, G. & Maeder, A. 2000, A&A, 361, 101 
Meynet, G. & Maeder, A. 2002, A&A, 390, 561 
Meynet, G. & Maeder, A. 2003, A&A, 404, 975 
Mihalas, D., Kunasz, P. B., & Hummer, D. G. 1976, ApJ, 206, 515 
Moffat, A. F. J., Schnurr, O., Chene, A.-N, et al. 2007, Highlights of Astronomy, 

14, 197 

Mokiem, M. R., de Koter, A., Evans, C. J., et al. 2007a, A&A, 465, 1003 
Mokiem, M. R., de Koter, A., Puis, J., et al. 2005, A&A, 441, 711 
Mokiem, M. R., de Koter, A., Vink, J. S., et al. 2007b, A&A, 473, 603 
Nahar, S. N. 1999, ApJS, 120, 131 
Nahar, S. N. & Pradhan, A. K. 1997, ApJS, 111, 339 

Najarro, F., Figer, D. E, Hillier, D. J., & Kudritzki, R. P. 2004, ApJ, 611, L105 
Nakamura, F. & Umemura, M. 2002, ApJ, 569, 549 
Niedzielski, A. & Skorzynski, W. 2002, acta, 52, 81 
Nugis, T. & Lamers, H. J. G. L. M. 2002, A&A, 389, 162 
Oskinova, L. M. 2005, MNRAS, 361, 679 

Oskinova, L. M., Hamann, W. R., & Feldmeier, A. 2007, ArXiv e-prints, 704 

Pauldrach, A. W. A., Hoffmann, T. L., & Lennon, M. 2001, A&A, 375, 161 

Petrovic, J., Pols, O., & Langer, N. 2006, A&A, 450, 219 

Pistinner, S. & Eichler, D. 1995, ApJ, 454, 404 

Puis, J., Springmann, U., & Lennon, M. 2000, A&AS, 141, 23 

Rauw, G., Vreux, J.-M., Gosset, E., et al. 1996, A&A, 306, 771 

Schmutz, W. 1997, A&A, 321, 268 

Schweickhardt, J., Schmutz, W., Stahl, O., Szeifert, T, & Wolf, B. 1999, A&A, 
347, 127 

Seaton, M. 1987, Journal of Physics B Atomic Molecular Physics, 20, 6363 
Seaton, M. J., Yan, Y., Mihalas, D., & Pradhan, A. K. 1994, MNRAS, 266, 805 
Seaton, M. J., Zeippen, C. J., Tully, J. A., et al. 1992, Rev. Mexicana Astron. 

Astrofis., 23, 19 
Spite, M., Cayrel, R., Plez, B., et al. 2005, A&A, 430, 655 
Springmann, U. 1994, A&A, 289, 505 

Stasihska, G., Grafener, G., Pefia, M., et al. 2004, A&A, 413, 329 

The Opacity Project Team. 1995, The Opacity Project, Vol. 1 (Institute of Physics 

Publications, Bristol, UK) 
Unsold, A. 1955, Physik der Sternatmospharen, 2nd edn. (Springer Verlag, 

Berlin) 

van der Hucht, K. A. 2001, New Astronomy Review, 45, 135 
Vink, J. S. & de Koter, A. 2002, A&A, 393, 543 
Vink, J. S. & de Koter, A. 2005, A&A, 442, 587 

Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2000, A&A, 362, 295 
Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2001, A&A, 369, 574 
Woosley, S. & Heger, A. 2006, ApJ, 637, 914 
Woosley, S. E., Langer, N., & Weaver, T. A. 1993, ApJ, 411, 823 
Yoon, S.-C. & Langer, N. 2005, A&A, 443, 643 



